function plotMeshWithCrack(model,elements,param)
mesh = model.getMesh();
ne = mesh.countFaces();

plotMesh(mesh,param);

for i=1:ne
    if elements(i).isCrack
        X1 = elements(i).X1;
        X2 = elements(i).X2;
        line([X1(1) X2(1)],[X1(2) X2(2)],'color','r');
    end
end

if cellfind(param,'positive domain')
    hold on
    plotPositiveDomain(model,elements);
end

end
